Trimers, molecules and polarons in imbalanced atomic Fermi gases 
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We consider the ground state of a single "spin-down" impurity atom interacting attractively with 
a "spin-up" atomic Fermi gas. By constructing variational wave functions for polarons, molecules 
and trimers, we perform a detailed study of the transitions between each of these dressed bound 
states as a function of mass ratio r — m^/mi and interaction strength. We find that the presence 
of a Fermi sea enhances the stability of the p-wave trimer, which can be viewed as a Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) molecule that has bound an additional majority atom. For sufficiently 
large r, we find that the transitions lie outside the region of phase separation in imbalanced Fermi 
gases and should thus be observable in experiment, unlike the well-studied equal-mass case. 



The spin-imbalanced Fermi gas has received much at- 
tention recently owing to its elegant realization in ultra- 
cold atomic gases. The ability to tune both the inter- 
species interaction and the spin polarization has allowed 
cold-atom experiments to access the rich phase diagram 
of the spin-imbalanced system However, as experi- 
ments become ever more precise, a challenge for theory 
is to go beyond conventional mean-field approaches and 
accurately determine the existence of exotic phases in the 
regime of strong correlations. 

One approach that has proven useful at high polariza- 
tions is to consider the problem of a single spin-down im- 
purity atom immersed in a Fermi gas of spin-up atoms 0- 
y. Such a scenario is just one example of the canonical 
"polaron" problem, the solution of which is used to con- 
struct the low-energy behaviour of the many-body sys- 
tem. Moreover, this limit of full polarization contains 
some of the critical points of the full zero-temperature 
phase diagram, e.g., it features the tricritical point that 
marks the existence of the spatially-homogeneous super- 
fluid phase for all polarizations [5|- Thus, an analysis 
of the high-polarization limit allows one to characterize 
parts of the topology of the whole phase diagram. 

An important feature of the single-impurity problem 
is that it can exhibit binding transitions where the im- 
purity changes its statistics and/or effective mass for a 
sufficiently strong attractive interaction. Thus far, the 
focus has been on equal masses mi — m-^, where it has 
been shown that the impurity undergoes a first-order 
quantum phase transition from a polaron (an impurity 
dressed with particle-hole excitations) to a molecule (an 
impurity bound to a single majority fermion) Here 
we demonstrate that a richer variety of phases can be ob- 
tained when the masses are unequal, with mi < m-f. In 
this case, the presence of a spin-up Fermi sea can favor a 
dressed molecule with nonzero ground-state momentum, 
corresponding to the FFLO superfluid phase in the limit 
of extreme imbalance. However, we find that the FFLO 
molecule generally prefers to bind another majority atom 
to form a p-wave trimer at zero momentum. Thus, as one 
increases the density of majority spins from the limit of 



zero density, one finds that the stability of the trimer is 
initially enhanced by the Fermi sea. 

In this Letter, we map out the ground-state phase dia- 
gram for the dressed impurity as a function of mass ratio 
r = m^/mi and interaction strength. We also examine 
the decay rates out of the various phases, after quenching 
across a phase transition. Lastly, we determine whether 
or not our dressed-impurity phase diagram is thermody- 
namically stable in the highly-imbalanced Fermi gas, by 
estimating the onset of phase separation in this limit us- 
ing the results of previous quantum Monte Carlo (QMC) 
simulations. In contrast to the equal-mass case, we find 
that the polaron-trimer binding transitions sit outside of 
the phase-separated region and should thus be experi- 
mentally accessible. 

In the following, we use the Hamiltonian for a two- 
component (we call the two species f, i) atomic Fermi 
gas interacting via a wide Feshbach resonance: 

=Xl'=kac],^Cka + |r 4^4,^Ck'+qiCk-qt , (1) 

kcr k,k',q 

where eko- = 2^ ^'^^ fi. = 1), is the system volume 
and g is the strength of the attractive contact interaction. 
The s-wave scattering length is then obtained via the 
prescription = — -t~ y^,t — r — i where the 'reduced 
mass' — = — — I — —, and A is a UV cutoff that can be 

nir m-^ ' 

sent to infinity at the end of the calculation. Note that 
since the J, impurity is distinguishable from the t Fermi 
sea, our results in the single impurity limit are relevant 
to both Fermi gases and Bose- Fermi mixtures. 

We construct variational wave functions for the single 
impurity atom immersed in a Fermi sea by considering 
different numbers of particle and particle-hole pair exci- 
tations upon the Fermi sea. Previous studies have shown 
this approach to be reasonably accurate Q. The nov- 
elties of our work are to thoroughly explore the case of 
unequal masses, to allow the dressed impurity to have 
non-zero momentum, and to consider the possibility of 
trimers as well as polarons and molecules. We also con- 
sidered tetramers at angular momentum L = 1, but they 
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did not bind anywhere for the parameters we considered. 

We adopt the following nomenclature for the different 
trial wave functions: a subscript n refers to a state with 
at most n operators acting on the non-interacting Fermi 
sea, and a bracketed momentum refers to the total mo- 
mentum of the state. For example, the polaronic state 
P3(Q) will contain the impurity atom, and at most one 
particle-hole pair on top of the Fermi sea, with a total 
momentum Q. In our explicit calculations we concen- 
trate on states with at most one particle-hole pair, since 
these are known to provide a good approximation for 
the impurity energy when m-j- — [6|t2| and going to 
higher order greatly increases the numerical effort. Thus, 
the wave function for the polaron is 0, 0] : 



|P3(Q)) = a^Q'ct 



\FS) 



(2) 



where \FS) is a Fermi sea of majority atoms, filled up 
to momentum kp^. The presence of a Fermi sea implies 
that the spin- up hole momentum = q < kp^ and the 
spin-up particle momentum k satisfies kp^ < k < K. 
Likewise, for the molecule, the wave function is 



k 

~ 4k'VQ+q-k-k'i4T4'T^<itl^^> • (3) 



k,k',q 

Finally, for the trimer, we set Q = and approximate the 
wave function as being insensitive to the hole momentum 
q by using its value at q = (this approach has been used 
for the polaron Q and molecule [1, Q , and can be shown 
to give an upper bound to the binding energy): 



^5(0)) = ^k,k.cU,_k.i4t4tl^^) + 

ki,k2 

??ki,k2,k,oc]j_k,_k,-k4,4it42t'=kt^qtl-F''5') • (4) 



E 

ki ,k2 ,k,q 



Here, we set the total angular momentum to be i = 1 
and choose = 0, in which case the bare part is 

Tkika = fci • zF{ki,k2, ki ■ k2) - h ■ zF{k2,ki,ki ■ k2) 

with any function F . Unlike the polaron and molecule, 
which both have L = in the ground state, we find that 
this L = 1 odd-parity trimer is always the lowest energy 
state throughout the portion of the {l/kp^as,r} phase 
diagram where the trimer is stable. We have not yet ex- 
plicitly looked at trimer states with nonzero momentum, 
but, as far as we are aware, there is no indication that 
the trimer has a lower energy at Q 7^ 0. 

In the limit kp^ — ^ 0, this L = 1 trimer becomes the 
3-particle bound state in a vacuum, for which analytical 
solutions have been found Here, it has been shown 
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FIG. 1: (Color online) The ground-state phase diagram as 
a function of mass ratio r and interaction strength l/kp-fas 
for the polaron (P3), molecule (M4) and trimer (Ts) wave 
functions. The FFLO region corresponds to M4 with non- 
zero momentum; the momentum of the FFLO molecule is 
approximately fc^t at the P3(0)-FFLO transition line and 
goes continuously to zero at the FFLO-M4(0) transition line, 
coming in as a square root near the transition. The T5-M4 
boundary (full [blue] line) approaches the 3-body transition 
rci — 8.17 in the limit l/kp^as — >■ 00, as expected. Above the 
dashed-dotted (r = rc2) line, the results for Ts become cutoff- 
dependent, and are therefore no longer universal. The shaded 
region marks where the system is unstable to phase separa- 
tion. The dashed (blue) line marks the transition line where 
a metastable polaron rapidly decays into a FFLO molecule 
which then decays into a trimer. 



that the trimer is bound relative to a molecule and an 
extra particle for r > rci = 8.17. However, for r > 
7'C2 — 13.6, the energy of the trimer is no longer finite in 
the limit A — >■ 00. This shows up as a wave function with 
weight at increasingly high momenta as one approaches 
rc2 from below. This critical rc2 is independent oi kp^, 
since it relies on high momenta. Thus, the results for the 
trimer are always cutoff dependent once r > rc2- 

We solve for the ground-state energy of a given normal- 
ized wave fimction by taking E = {ip\H ji/') and then 
minimizing the energy E with respect to the amplitudes 
a, (3, 7, S, T and/or 77. One can generically sum over 
one of the momenta, and solve a self-consistent equa- 
tion for the resulting function, where symmetries play 
a crucial role in simplifying the calculation. For exam- 
ple, for the trimer with L = 1 and = 0, we define 
V Ska ^k,k2 = G{k)k-z, and then derive a self-consistent 
equation for G. 

By determining the minimum energies of the wave 
functions /3(Q), A/4(Q), 15(0), we constructed an ap- 
proximation to the ground-state phase diagram of the 
impurity atom, as depicted in Fig. [TJ Here, we find that 
the polaron phase always has its lowest energy at zero 
momentum. By contrast, a very small region in Fig. [1] 
exists where the molecule has nonzero momentum in the 
ground state (a FFLO molecule). 
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The FFLO phase has a very natural interpretation in 
this hmit: having a molecule with Q = requires the im- 
purity atom to have momentum k > kp^, which is kinet- 
ically disfavored as r increases. Thus, a FFLO molecule 
appears once the impurity's momentum drops below kp'f- 
However, as r is increased farther, it quickly becomes fa- 
vorable for the FFLO molecule to bind yet another ma- 
jority fermion and form a p-wave trimer at zero total 
momentum. Indeed, at large enough l/^kpfUs), the tran- 
sition is directly from the Q = molecule to the trimer, 
without any intervening FFLO phase. A surprising re- 
sult is that as one approaches unitarity, the Fermi surface 
actually favors the trimer phase over the molecule phase 
and we have a direct transition from polaron to trimer. 
We expect this to be a robust result since our approxi- 
mation for T5 underestimates the trimer binding energy. 

How are these phases distinct from one another? In the 
polaron phase, the ground state wavefunction (including 
dressing with an arbitrarily large number of particle-hole 
pairs) has a nonzero weight Z = |q;^'*^-'P for the "bare" 
impurity atom. This weight appears to remain nonzero in 
the thermodynamic limit only at momentum Q = 0: the 
effective mass m* of the polaron appears to be positive 
whenever the polaron is the ground state in our phase 
diagram, thus a polaron with a nonzero Q has higher en- 
ergy by / (2m*) and is unstable to emitting a particle- 
hole pair and scattering to lower momentum. We may 
consider this weight Z to be the "order parameter" for 
the polaron phase. 

In the molecule phase, we have Z = 0, but a fully- 
dressed ground state will have a nonzero weight Zm — 
X^k l7k'^''l^ for the bare molecule at one or more values 
of Q. The FFLO molecule phase is where this nonzero 
weight appears at Q > 0. Thus again we can consider this 
nonzero Zm, as well as the momentum Q where it occurs, 
as the "order parameters" of the molecule phases. Sim- 
ilarly, the trimer phase is where the ground state has a 
nonzero weight Zx = X^k k I''kik2p for the bare trimer. 

Next, we examine the nature of the binding transitions 
for the single impurity. For the case of equal masses, it 
is known that the matrix elements connecting the po- 
laron and molecule go to zero at the transition point due 
to the transition involving a "decay" into at least four 
particles Q. Thus, one can have a long-lived molecule 
(polaron) existing on the polaron (molecule) side of the 
transition, with the lifetime F"^ ^ {AE)~^^'^ diverging 
faster than the inverse of the energy difference AE as 
one approaches the unbinding transition [lH . This first- 
order binding transition in the presence of the Fermi sea 
occurs when the momenta of the two phases do not differ 
by A;Ft- For example, when the molecule and polaron 
both have their minimum energy at Q = 0, the molecule 
at Q = cannot simply unbind into a polaron at Q = 
plus a low energy spin up particle at the Fermi momen- 
tum kp-f, without violating momentum conversation. In- 
stead, it must create a particlc-holc pair, thus leading to 




FIG. 2: (Color online) Schematics of two different scenarios 
for a molecule unbinding into a polaron -I- particle. The solid 
(red) lines represent the molecule dispersion E{Q) and the 
shaded regions correspond to the polaron -I- particle (two- 
body) continuum. When both the molecule and polaron have 
their minimum energies at Q = (left), the transition is 
first-order. However, we have a continuous transition (where 
the bound state fully "mixes" with the continuum) when the 
molecule has ground-state momentum Q — kpf (right). 



the 4-particle process mentioned earlier and a first-order 
transition. Referring to Fig. [21 this follows from the dis- 
placement of the molecule dispersion from the polaron -|- 
particle continuum by kp-f. The same considerations ap- 
ply to the zero-momentum trimer unbinding into a zero- 
momentum molecule plus a particle. 

However, if the molecule has its energy minimum at 
Q = kp^, then it can simply unbind in a continuous fash- 
ion into either a zero-momentum polaron plus a particle 
(sec Fig. [2]), or a zero- momentum trimer plus a hole. This 
behaviour also shows up in the decay rates of the excited 
states on either side of the transition, which we can esti- 
mate using a simple Fermi's golden rule like in Ref. 
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We find that the unbinding of a FFLO molecule into a 
trimer or polaron has decay rate F ^ AE as AE 0, 
where, once again, AE is the energy difference between 
the two phases near the transition. Therefore, this is a 
marginal case where the lifetime F"^ diverges as fast as 
the inverse of AE. The decay of a trimer or polaron into 
a FFLO molecule is even faster, with F being finite as 
AE — )■ 0, since the momentum of the final-state FFLO 
molecule can lie anywhere on the Fermi surface and thus 
the phase space is enlarged. Finally, the transition from 
the zero-momentum molecule to the FFLO molecule is 
not an unbinding transition and we find that Q generally 
moves continuously away from zero at this transition. 

A direct transition from the trimer to the polaron and 
vice versa would be first-order, since this is a three-body 
decay: the trimer (polaron) "shedding" two particles 
(holes). In this case, the decay rate is always F ~ (AE)'^ 
close to the transition, similar to the behavior of a Fermi 
liquid quasiparticle. This implies that a quench at fixed r 
from the polaron to trimer phase can lead to a metastable 
polaron that can become unstable to forming a FFLO 
molecule, as shown in Fig. [T] If the difference between 
the FFLO and trimer energies is small here, then there 
may exist a metastable FFLO phase, a scenario which is 
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worthy of further investigation. 

In experiments on highly-imbalanced Fermi gases, one 
has a nonzero number density of minority atoms, and 
the system must be stable against phase separation in or- 
der for the single-impurity transitions to be observable. 
For equal masses, QMC predicts that the single -imp urity 
transition occurs in the phase-separated regime [l2| , and, 
indeed, the experimentally measured disappearance of 
the polaron [13[ agrees well with the onset of phase sep- 
aration between superfluid (SF) and normal (N) phases, 
rather than with the polaron-molccule transition for a 
single impurity. To estimate the onset of phase sepa- 
ration for general mass ratios, we impose the following 
coexistence conditions on the pressures P and chemical 



potentials fi in each phase: Psf ~ Pn, IJ'a — t^a ■ ^.nd 
/ij^ = El, , where Ei, is the binding energy of the impurity 
(a polaron, molecule or trimer) immersed in the Fermi 
gas. Fortunately, for mass ratios r > 1, the superfluid 
is likely to be unpolarized (n-f = n^) or only weakly po- 
larized at this onset point, and thus we can exploit the 
QMC equation of state for the unpolarized superfluid. 
At unitarity, this equation of state has been shown to 
be relatively insensitive to mass ratio while in the 
BCS and BEC limits, the equation of state is that of a 
weakly-attractive Fermi gas and a weakly-repulsive Bose 
gas, respectively. Thus, we have pressure and average 
chemical potential: 



PSF = 



3/2 



_5/2 



157r2 



^J■SF 



g{\/kFas,r) (5) 
Ef f{l/kFas,r) (6) 



where the functions g{l/kpas, r) and /{l/kpUs, ?') arc de- 
termined via numerical interpolation between the known 
limits, while kp = (GTr^n^^)^/^ and ep = kl^/2'mr. 

Clearly, the pressure and spin-up chemical potential in 
the fully-polarized normal phase are known exactly: 



Pn - 
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-5/2 
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(7) 
(8) 



where ep^ = {kp^)"^ /2mr. To determine /ij^, we estimate 
the binding energy Ei, using the P3 wave function. Note 
that at large r, where the polaron is no longer the ground 
state, this will be an underestimate, and so the calcu- 
lated onset position 1/kp^as will be lower than the exact 
result in this limit. By applying coexistence conditions, 
we arrive at the set of equations 



kp 



(1 + r) g{l/kpa,,r) 



(9) 



5/2 
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iyk,a..r)^il,rf'^[mi^^ (10) 



For equal masses, our calculation gives l/kp^ag ~ 0.75, 
which agrees well with fixed- node QMC calculations [l^l ■ 

This onset line determines the region of phase separa- 
tion near unitarity and this is plotted in Fig. [T] For suffi- 
ciently large r, wc see that the polaron-trimer transition 
extends well outside of the regime of phase separation, 
and thus should be observable in the polarized gas. This 
allows one to investigate first-order binding transitions 
and the possible metastable states discussed earlier. 

A remaining question is what happens to the trimer 
phase when there is a finite density of spin-down atoms. 
In one dimension, it is known that a trimer phase exists 
for n \ /n-^ = 1/2, provided the interactions are sufficiently 
large [l5|. Here, wc expect a Fermi liquid of trimers for 
low densities n^/n^ <^ 1/2. However, it is possible that 
the FFLO phase will eventually win over the trimer phase 
as ?^4./^^.-|■ is increased. One may also have a mixture of 
trimers, polarons and/or molecules as we approach the 
single-impurity binding transition. 

Experimentally, one can explore our phase diagram us- 
ing two different atomic species with unequal masses or 
by artificially increasing the effective mass of the spin- 
up atoms with a spin-dependent optical lattice. "'"K- 
^Li mixtures are the favored choice, and, in Fig. [U their 
mass ratio (r ~ 6.7) just touches the bottom parts of the 
trimer and FFLO regions. Since each phase is charac- 
terized by a different effective mass, one way of distin- 
guishing these phases is to determine the effective mass 
using the low-lying compression modes of the gas, as in 
Ref. . Another possib ility is to directly measure the 
"order parameter": Ref. [l3| has successfully measured 



Z in the polaron phase using RF spectroscopy, but it re- 
mains an open question whether or not this method can 
be extended to study Zm or Zt- 
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